New possibility of the ground state of quarter-filled one-dimensional strongly 
correlated electronic system interacting with localized spins 
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We study numerically the ground state properties of the one- dimensional quarter- filled strongly 
correlated electronic system interacting antiferromagnetically with localized 5* = 1/2 spins. It is 
shown that the charge-ordered state is significantly stabilized by the introduction of relatively small 
coupling with the localized spins. When the coupling becomes large the spin and charge degrees of 
freedom behave quite independently and the ferromagnetism is realized. Moreover, the coexistence 
of ferromagnetism with charge order is seen under strong electronic interaction. Our results suggest 
that such charge order can be easily controlled by the magnetic field, which possibly give rise to the 
giant negative magnetoresistance, and its relation to phthalocyanine compounds is discussed. 

PACS numbers: 71.20.Rv, 71.30.+h, 74.70. Kn 



In the quest of conducting states in molecular solids 
for the past thirty years, charge transfer complexes with 
the particular 2:1 composition of two different kinds of 1~L 
molecules, usually donor and acceptor, have played cru- 
cial roles, e.g. TMTSF2X is the classical example to 
show the first organic super conduct ivity(SC) and very 
rich phases on the plane of pressure and temperature jl|. 
The donor molecules have an average valence of +1/2 
forming a quarter-filled tt electronic system. There, an 
interesting competition and coexistence in the types of 
insulating states are seen under the strong electronic in- 
teraction depending on the degree of dimerization; charge 
ordering(CO) and dimer-Mott expected in the absence of 
(or under weak) dimerization and at under strong dimer- 
ization, respectively 0. Rich physics in these classes 
of materials have now been explored. In addition, an- 
other family including both quarter-filled 7r-electrons 
and localized spins showed up. For example, in two- 
dimensional A-BETSoFeCUj^ and one-dimensional(lD) 
EDT-TTFV02FeBr4 3 each acceptor molecule has 
6'=5/2(Fe^+), the same population number as elec- 
trons, whereas in TPP[Fe(Pc)(CN)2] 2 (phthalocyanine 
compound) each donor molecule has 5'=l/2(Fe^+) 5], 
twice the electron number. The combination of localized 
spins with the strongly interacting quarter-filled electrons 
offer even more interesting possibilities like field-induced 
SC[3|, ferrimagnetismQ, etc. 



In this paper we focus on the ID TPP[Fe(Pc)(CN)2]2 
and study a simplest model (extended Hubbard-Kondo 
model) including an antiferromagnetic Kondo coupling, 
J, between the conduction electrons and the localized 



spins, whose Hamiltonian is given by, 

{'ij)cr (ij) 3 j 

Here, Cjcr, and Sj are annihilation, number and spin 
operators at site j for conduction electrons, respectively. 
Sj represent spin operators of the localized spins with S = 
1/2, and {ij) denote the nearest-neighbor (nn) pair sites. 
We do not take account of the effect of dimerization, 
which holds in TPP[Fe(Pc)(CN)2]2. We focus on the 
quarter-filling of the conduction electrons where many 
organic conductors as well as the topic material reside. 

For the case with J = 0, it has been confirmed that the 
CO insulating (C 01) state is realized at around U > At 
and V >2ty^ 8]. On the other hand, usual Kondo lattice 
model does not include the Coulomb interactions, i.e., 
U = V = {). In this case, when the Kondo coupling is larger 
than some critical value, Jc{^)^ dependent of the electron 
filling, J 9, 10], a ferrimagnetic metallic state is realized, 
i.e., all the conduction electrons form singlets with the 
localized spins and the remaining localized spins align in 
one direction. We call this state as ferromagnetic in the 
following. This model has already provided some pictures 
of the interplay of charge and spin degrees of freedom for 
the heavy fermion systems [Of. Still, the strong electronic 
interaction adopted in the present model would make the 
physics richer. In this context, we will study the interplay 
between the COI and ferromagnetic metal(FM). 

First we examine whether the COI state of the ex- 
tended Hubbard model (J = 0) survives in the fi- 
nite J region. This can be examined by calculat- 
ing the charge gap, Ac(iV) = EN^Ne ^2) ^ EN^Ne - 
2) — 2EN{Ne)^ and extrapolating it to infinite system 
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FIG. 1: (a) Charge gap Ac as a function of at fixed value 
of J/t — 1 and U/t = 8 for several choices of V/t. Typical 
truncation error is 10~^-10~^ with a system size up to A/" = 
200. Symbols at 1/N = are the extrapolated values, Ac(oo). 
(b) Ground state phase diagram on the plane of U/t and V/t. 
The region above the solid (dashed) lines for J/t= 1 and 5 
{J/t — 0, oo) are the COI state where Ac(cxo) remains finite. 

size. Here, EN{Ne) is the ground state energy of 
system size N with the electron number Ne. Fig- 
ure nja) shows some examples of the A^-dependence 
of Ac obtained by the density matrix renormalization 
group (DMRG) method |12i|. We adopt a so-called non- 
Abelian algorithm [15 \l4, Il5j with open boundary con- 
dition. In this algorithm, we deal with only four states 
per single site, (5, ne) = (l, 1), (1/2, 0), (1/2, 2), (0, 1), clas- 
sified by the total spin, s, and the electron number, ng, 
e.g. (1,1) represents a triplet state formed by a localized 
spin and a conduction electron. The use of this reduced 
basis provides a dramatic performance of calculation as 
well as perfect conservation of the SU(2) symmetry which 
is quite easily broken in the conventional DMRG. 

From Fig^a), the charge gap at ^00 is estimated 
by fitting the data to Ac{N) = Ac(oo) + A/iV + B/N^, 
where A and B are constants. The metal-insulator (MI) 
phase diagram at several fixed values of J/ 1 derived from 
this extrapolation is shown in Fig. 1(b). We use the 
lowest-energy state irrespective of the total spin. For 
J/t = 0, Fig. 1(b) reproduces the previous results by the 
exact diagonalization in small clusters 0, and the insulat- 
ing state is the CO state. When J/t = 1, this COI phase 
is significantly stabilized. This is rather surprising since 
the charge degrees of freedom is easily controlled by a 
small amount of interaction with localized spins. By fur- 
ther increasing J/t, the phase boundary asymptotically 
converges to the V/t = l line in the limit of J/t = 00. 

There are several physical reasons of this enhancement 
of COI state by the Kondo coupling. Basically J intro- 
duces a tendency to form a singlet. (I) Firstly, this re- 
duces the effective hopping of the conduction electrons. 
For example, at J/t 1, hopping occurs only when one 
singlet is broken at a site and then a new singlet is formed 
on its nn site. In this case, the hopping matrix element is 



teff = t X (1/V2)^ = tf^j and thus V/t^f^ increases from 
the original value of V/t^ which favors the CO state. (II) 
Another reason is that the singlet formation disfavors 
the double occupancy of the conduction electrons. This 
works cooperatively with U. (Ill) Finally, a singlet state 
can gain a self-energy by forming a virtual triplet state 
on the nn sites. This occurs only when the nn site is not 
occupied by another conduction electron. This plays a 
similar role as V. All these effects extend the CO region 
towards the smaller U/t and V/t region. 

These phenomena become clearer if we study the 
strong-coupling region of J/t ^ 1 and U/t^ 1. In this 
region, each single electron forms a singlet with a local- 
ized spin on the same site. We can regard this singlet as 
a vacant site and the remaining unpaired localized spin 
as "particle" which cannot have double occupancy. Such 
"particle" is denoted by dja- in the following. The hop- 
ping of conduction electrons or singlets is interpreted as 
that of the "particles" in the inverse direction. For small 
t/J and t/[/, we take account of the second order per- 
turbation, and the effective Hamiltonian becomes, 

JO- 

+jjd;+2adjanj+i j + h.c. 

jrr' 
3 

where and denote the spin and number operator 
of the "particles" , respectively. The second and the third 
terms in the first parenthesis give modification of the hop- 
ping integral, and the next term with t^/ J causes the spin 
interaction which will be discussed shortly. The last term 
represents the modification of the V term as mentioned in 
the reason (III). In the limit of J = (X), [/ = oo, TYeff is re- 
duced to the U = 00 limit of the extended Hubbard model 
with teff = because the double occupancy of "parti- 
cles" is not allowed. Since the nn interaction between 
"particles" is still equal to the MI phase boundary is 
^/^eff = 2, i.e., ]//t = l, consistent with Fig.QJb). 

Let us switch to the magnetic properties. Figure 2(a) 
shows the ground state phase diagram classified by the 
total spin, obtained in our DMRG. The phase transi- 
tion from a paramagnetic to a ferromagnetic state takes 
place at around J/t ~ 1 — 2 depending on the values of 
U/t and V/t. At[/ = ]/ = 0, a phase transition takes 
place at Jc/t ~ 1.6, in concurrence with the results of 
R.ef . [loj . The phase transition is the second order in the 
sense that the total S changes successively from to N/ A. 
However, the region of intermediate S is small and de- 
pends quite much on the numerical condition as in the 
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FIG. 2: Ground- state phase diagrams in the space of (a) U, V, 
and J, and its cross-section planes {h)U-V for several values 
of J/t, {c)J-U for several values of V/t, and (d)!/- J for several 
values of U/t. Typical truncation error is 10~^. 

case with U = F = O0, so that we determine the phase 
boundaries from the energy difference between S = and 
S = N/A states. We also confirm that the spin gap is 
absent within the numerical error at the representative 
parameter point of each phases. Figure Efb) shows the 
ground state phase diagram on the plane of U/t and V/ 1 
which is to be compared with Fig.^b). Apparently, the 
magnetic and MI phase boundaries are not associated 
with each other and the spin- and charge-degrees of free- 
dom behave independently. Although a charge gap opens 
above some Vc, this does not affect the spin degrees of 
freedom much because the charge gap opens quite slowly. 

The detailed cross-section phase diagram on the plane 
of U/t and J/t is given in Fig.E^c). The phase boundary 
moves towards smaller J/t when finite U/t is introduced. 
This is because the singly occupied states of conduction 
electrons are stabilized hy U/t. From Figs.El^c) andElJd), 
we can see that the introduction of V/t slightly stabilizes 
the ferromagnetic phase when U/t is large, whereas the 
opposite tendency is observed when U/t is small. When 
V/t is larger than /7/t, two electrons occupy a single site 
rather than to stay apart next to each other, and an in- 
stability towards superconductivity takes place 0. Such 
double occupancy hinders the ferromagnetism. This is, 
however, quite an unrealistic situation. 

In order to study the ferromagnetism in the large U/t 
region, it is convenient to use the effective Hamiltonian 
of eq.(|2|). As in the large-[/ Hubbard model, the ground 
state wavefunction can be expressed as. 



where a spinless fermion wavefunction, '^srC^eff, Kff), 
represents the charge degrees of freedom and the 
spin degrees of freedom on a squeezed chain 16]. In the 
limit of J = oo, = 00, the spin degrees of freedom 
are degenerate as in the U = 00 Hubbard model. Then 
the second-order terms in eq.© can be considered as 
perturbations 17]. The first parenthesis and the last term 
in eq.(|21) modifies t and V of the spinless fermions into 
teff and , respectively. The degeneracy of the spin de- 
grees of freedom is lifted by the remaining terms. These 
terms give an effective ferromagnetic Heisenberg spin in- 
teraction, Jeff • 'S'i+i, in the squeezed spin chain with 



Jeff = y(^./+2^j^i+i)sF < 0, 



^ =^SF(teff, V;ff)<l>H, 



(3) 



where (• • •)sf represents the expectation value in '^sfE^- 
This leads to the ferromagnetic ground state of <I>h. 

Let us discuss the V-dependence of Jeff. When V or 
Feff is introduced it reduces the nn population of "par- 
ticles", which suppresses Jeff. When the CO becomes 
dominant, we expect a very small Jeff. However, in- 
finitesimal Jeff always stabilizes the ferromagnetic state 
because of the degenerate perturbation. Actually, we 
find, Jeff = -Jt^/2V'^{V + J)^ in the large V-region by 
considering the fourth order hopping process of Fig.IHlJa); 
we assume the two localized spins with a singlet in be- 
tween. The electron hops to its nearest neighbor to form 
a singlet or a triplet as a virtual state. The triplet forma- 
tion in the middle of the process (after hopping back to 
the original site) allows the mixing between left and right 
neighboring localized spins to give Jeff. Among the four 
different paths in Fig.[2fa), the path including singlts in 
the first and third virtual state corresponds to the process 
included in eq. (4). Therefore, the ferromagnetism due 
to Jeff in eq. (4) naturally continues to that of the large 
V-region. We calculate the nn and next-nearest neighbor 
(nnn) correlation functions between charges and localized 
spins as shown in Fig.|3Jb); nn values are suppressed by 
V in contrast to the nnn ones, which proves the gradual 
crossover from the ferromagnetism by Jeff to that by Jeff. 

The ground state phase diagram of eq.Q is summa- 
rized in Fig.^ at small J/t the phase diagram consists of 
paramagnetic metal(PM) and paramagnetic COI(PCOI) 
at small and large U/t and V/t^ respectively. With in- 
creasing J/t, PCOI phase expands. At the same time, 
the ferromagnetic phases (FM and FCOI) descend sig- 
nificantly from large /7/t-part at around J/t^l — 2. 

Let us discuss the possible implications of the present 
theoretical results to experiments in TPP[M(Pc)(CN)2]2, 
M=Fe,Co. In the case of Co-salt which does not have lo- 
calized spins, i.e. J=0, a weak sign of CO is observed 
in the NQR and the resisitivity shows semiconducting 
temperature dependence with a small activation energy 
(Aft ~ 10~^eV). This indicates that Co-salts are in the 
CO region but close to the MI phase boundary as shown 
in FigElin view of U/t -6-8 and V/t -3 with t -O.leV 
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FIG. 3: (a) Fourth order perturbation process at J,U,V ^ t 
in eq.(QJ. Filled and open big circles denote the singlet and the 
triplet, respectively, and the site without them is composed 
of a free localized spin without electrons, (b) Localized spin- 
spin and charge- charge correlation functions of nn and nnn 
sites averaged over the whole system as a function of V/t at 
the fixed value of U/t = 8.0 and J/t = 1.0. 




FIG. 4: Summarized phase diagram on the plane of U/t 
and V/t. Black circle represents the possible location of 
TPP [M(Pc) (CN)2] 2 (M=Co,Fe) . 



to make Fe-salt a POOL Such small J/t does not lead 
to ferromagnetism which also agrees well with the exper- 
iment. We speculate that the weak ferromagnetism at 
T ~5K is due to the canted magnetic moments induced 
by impurities in the antiferromagnetic background. 
Interestingly, in the Fe-salt, a giant negative magnetore- 
sistance(GNMR) is reported below T-SOkTs]. This will 
simply be understood in the present framework(eq.(P)) as 
due to the suppression of CO phase. When the magnetic 
field is applied, the total Sz increases and each site will 
have the larger population of localized 52=1/2 (with- 
out electrons) as well as of Sz=l triplets composed of 
one eletron and a localized spin, whereas the number of 
singlets with Sz=0 decreases. Since the electron hop- 
ping from 6';2=l-site to Sz=l/2-site gives teff=l, which 
is larger than those assosiated with Sz=0-site{tef[ = t/^/2 
or t/2), y/teff is reduced and the charge gap will be sup- 
pressed. This process of GNMR is quite different from 
the conventional double-exchange mechanism 19]. 
In conclusion, the CO in the quarter-filled strongly cor- 
related electronic system becomes stable when combined 
with the localized spins, and this CO could be easily con- 
troled by the magnetic field. 



given by the extended Hiickel calculation. The mag- 
netic susceptibility behaves as that of the ID Heisen- 
berg 5'=l/2 spin system. Here, the antiferromagnetic 
coupling is expected to be several factors larger than 
the one between spins of completely localized charges 
(Jaf = t^/UV'^ ~ lOK), because of a certain degree of 
delocalization in the proximity of the MI boundary. 
As for the Fe-salt with localized 5'=l/2 spins, the re- 
sistivity is semiconducting with larg er /A(2 IQ-^eV, 
which indicates the stronger CO in Fe salt than in Co- 
salt. This is consistent with the present results that 
the coupling to the localized spins generally enhances 
the charge gap. The magnetic susceptibility in Fe-salt 
shows large anisotropy between those parallel(x||) and 
perpendicular (x^) to the ID chain; X|| behaves similar 
to the Co-salt case, while x± obeys the Curie-law at 
T > T/v ^25K and suddenly saturates to become weakly 
T-dependent at T ~ 5K-25K, followed by a weak ferro- 
magnetism at T <5K. Noting that magnetic easy axis of 
Fe ions is esentilally perpendicular to the ID chain, we 
may assume that X|| and x± are originated mainly from 
the electrons and the Fe spins, respectively. Hence, the 
similarity of X|| with that of the Co salt is understood, 
while the behavior of x± will be interpreted as follows; 
the saturation below Tn suggests that the antiferromag- 
netic correlation developes below T/v, which is due to the 
onset of TT-d interaction roughly estimated as J ~ T^. In 
our model, this J/t ~0.025 actually pushes enough the 
MI phase boundary to the region of smaller U/t and V/ 1 



[1] D. Jerome, Science 252, 1509, (2004). 

[2] H. Seo, C. Hotta, H. Fukuyama, Chem. Rev. 104, 5005, 

(2004), and the references therein. 
[3] H. Kobayashi, A. Kobayashi, T. Cassoux, Chem. Soc. 

Rev., 29, 325, (2000). 
[4] S. Noguchi, A. Matsumoto, T. Matsumoto, T. Sugimoto, 

T. Ishida, Physica B 346-347 397 (2004). 
[5] T. Inabe, H. Tajima, Chem. Rev. 104, 5503 (2004). 
[6] S. Uji, H. Shinagawa, C. Terakura, T. Terashima, T. Yak- 

abe, Y. Terai, M. Tokumoto, A. Kobayashi, H. Tanaka, 

and H. Kobayashi, Nature 410, 908 (2001). 
[7] F. Mila and X. Zotos, Europhys. Lett. 24, 133 (1993). 
[8] H. Yoshioka, M. Tsuchiizu and Y. Suzumura, J. Phys. 

Soc. Jpn. 69, 651 (2000). 
[9] H. Tsunetsugu, M. Sigrist, and K. Ueda, Phys. Rev. B 

47, 8345 (1993). 
[10] I. P. McCulloch, A. Juozapavicius, A. Rosengren, M. Gu- 

lacsi, Phys. Rev. B 65, 052410 (2002). 
[11] G. Aepph, Z. Fisk, Comments Cond. Mat. Phys. 16, 155 

(1992). 

[12] S.R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[13] I. McCulloch, M. Gulacsi, Europhys. Lett. 57, 852 
(2002). 

[14] In this algorithm at S t^O, the newly obtained basis after 
the truncation is a linear combination of the old basis 
with several total spins, i.e. the new density matrix has 
off-block-diagonal elements between different total spins. 
However, neglecting the off-block-diagonal elements does 
not alter the results within the numerical error 13]. 

[15] The number of states at each block, m, are kept up to 
200 which is enough to suppress the maximum truncation 
error to less than the order oflO~^-10~^ in the worst case. 

[16] M. Ogata, H. Shiba, Phys. Rev. B41, 2326 (1990). 



5 



[17] H. Shiba, M. Ogata, Int. J. Mod. Phys. B5, 31 (1991). (1999); K. Kubo, N. Ohata, J. Phys. Soc. Jpn. 33 21, 

[18] N. Hanasaki, H. Tajima, M. Matsuda, T. Naito, T. Inabe, (1972). 

Phys. Rev. B 62, 5839 (2000). 
[19] A. Moreo, S. Yunoki, E. Dagotto, Science 283 2034 



